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Abstract 

In this letter we calculate the exact partition function for free bosons on the plane with lacunae. 
First the partition function for a plane with two spherical holes is calculated by matching exactly 
for the infinite set of Wilson coefficients in an effective world line theory and then performing 
the ensuing Gaussian integration. The partition is then re-calculated using conformal field theory 
techniques, and the equality of the two results is made manifest. It is then demonstrated that 
there is an exact correspondence between the Wilson coefficients (susceptabilities) in the effective 
field theory and the weights of the individual excitations of the closed string coherent state on the 
boundary. We calculate the partition function for the case of three holes where CFT techniques 
necessitate a closed form for the map from the corresponding closed string pants diagrams. Finally, 
it is shown that the Wilson coefficients for the case of quartic and higher order kernels, where 
standard CFT techniques are no longer applicable, can also be completely determined. These 
techniques can also be applied to the case of non-trivial central charges. 



I. INTRODUCTION 



The partition functions for conformal field theories on surfaces with boundaries are of 
interest in diverse areas of physics. In soft condensed matter physics it is of interest to 
calculate partition functions on the Euclidean plane with compact boundaries for the purpose 
of understanding fluctuation induced forces |2j. In string theories, amplitudes correspond 
to such partition functions and one must sum over all moduli and Euler characteristics. 
However, to calculate a particular gauge slice is not always simple. Powerful methodologies 
in boundary conformal field theory (BCFT) jl| have been utilized to determine the spectrum 
of primary operators in Gaussian as well as non-Gaussian theories. These methods utilize 
the state operator correspondence to implement the boundary conditions via the insertion 
of boundary operators. 

In t 
(EFT) 

case of partition functions on the Euclidean plane with compact boundaries. The EFT 
approach is similar to BCFT in that the effects of the boundaries are accounted for via 
operator insertions. However, we will make no use of the state operators correspondence, 
and as such, the technique will not be limited to conformal field theories. 
We will begin by considering the free bosonic action 

S* = \ f d 2 r 3^ , (1) 

on a plane with holes which satisfy Neumann boundary conditions. This problem is identical 
to finding the partition function for a film embedded with rigid disks which are allowed to 
bob but not tiltQ, where the notion of tilting is given meaning by a particular choice of 
gauge. In particular the scalar field partition function will correspond to the Monge gauge 



lis paper we will take a different approach based on world-line effective field theories 
3] as applied to statistical systems {4, 5]. In particular, we will be interested in the 



choice jsl where the height of the field is measured with respect to a Euclidean flat base 
plane. In the film model there are non-linear corrections which are there as a consequence of 
reparameterization invariance, but it can be shown that these corrections are suppressed 4^ 
by powers of k/(3 where k is the modulus characterizing the surface energy. The fluctuation 
induced forces in these systems are relevant to soft condensed matter problems 



1 This is not true for all observables but will hold for the partition source-free partition function. 
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In |4( it was shown that it is possible to utilize the world-line effective field theory de- 
veloped for studying black hole inspirals js| to study the Casimir forces on two dimensional 
surfaces (the dimensionality is actually irrelevant). The idea is to work in the point particle 
approximation, and then account for the finite size (read boundary conditions) by adding 
an infinite tower of higher dimensional operators to the action. For defects with only lin- 
ear responses to external perturbations these terms are all quadratic in the field, and their 
dimensionality is fixed by the number of derivatives. Power counting shows that non-linear 
interactions on the world-line generate corrections which are suppressed at high temper- 
atures. In the limit where the size of the objects (R) is small compared to the distance 
between them r, truncating the infinite sum induces errors of order (R/r) T , where r is the 
the number of derivatives in the first term dropped in the series. The coefficients of these 
operators (Wilson coefficients) correspond to generalized polarizabilities and are fixed by 
performing a matching calculation between the full theory and the effective theory. The full 
set of Wilson coefficients can be calculated by first matching the response of the individ- 
ual holes to multipole fields in the full and effective theory. We will then show that these 
matching coefficients can equally well be extracted by determining the boundary operators 
in the Hamiltonian version of the theory. Comparing the two procedures illuminates a map 
between closed string excitations and effective polarizabilities. The model is then solved by 
summing the series for the higher dimensional operators into a non-local operator which can 
be used to write down a closed form expression for the partition function. 

In the effective theory the action for for n-bodies may be written as 

S = S + i^QC^....^)^....^), (2) 

i 

where the second term is localized at the origin. Operators involving traces over indices 
are redundant in this theory as they can be traded for other operators via a canonical 
transformation. For non-spherically symmetric objects the polarizibilities will transform 
non-trivially under the two dimensional rotation group, but are singlets for the spherical 
case of interest here. These coefficients are fixed by a simple matching procedure whereby we 
place the objects in a background field and calculate its response. For a given polarizability 
of a given order we may extract the corresponding Cj by placing the object in a multipole 
field of the same order. 



3 



To utilize holomorphy we match in the complex plane a general n-pole background field 
takes the form 

<p b9 (z,z) = ^(z n + z n )- (3) 
We will impose Neumann boundary conditions which preserve the conformal symmetry. 

d r <f>{z) \ R = (4) 

This choice means that the first non-zero (7, will be C*2, a dipole polarizability. In the full 
theory we solve the rudimentary boundary value problem. The induced field is given by 

^ = ^"(1 + 1). (5) 

Expanding around the background field ([3]) utilizing the Greens function in complex 
coordinates 

G(M) = ^0og(*)+log(*)) (6) 
and calculating the one point function leads to the EFT result 

We may then extract the matching coefficient to be 

° n = ~n\(n-\)\2^ R2n - (8) 

This form of matching coefficient can also be obtained from perspective of closed string 
theory. If we conformally map an isolated hole into an infinite cylinder, and consider the 
direction along the length as being time-like, then the partition function corresponds to a 
closed string propagator emanating from a boundary state Q. The boundary states are 
annihilated by the constraint 

<9i0 | V) = (9) 
where the world-sheet is parameterized by (£71,02) an d o\ is the time like direction. 

Our gauge choice (Monge) corresponds to choosing mode oscillators in the string expan- 



sion to be along one space like direction in the embedding space. In terms of the mode 
operators the constraint [3] reads 



(a n + a- n ) | V) = (10) 
where the (standard) normalization is fixed via 

/ ,\ 1/2 oo 



2 / ^ z r 

7 m=— oo 

[a m , a n ] = [a m , a n ] = m5 m - n , (12) 

and 

The solution to this constraint is 

| V) oc e " E ^i ( p a - p5 - p) | 0). (14) 



Using the fact that 



2 \ 1/2 i 



a -™=l,^J i^Ty 9 "^ < 15 > 



the state is written as 



| V) oc e~ E ~ =1 ^(l^^^) | 0). (16) 
rescaling the fields so that they are normalized according to (TjQ) leads to 



ip) oc Exp 



0). (17) 



Taking this as an operator insertion in the path integral reproduces the Wilson coefficients 
C n . It is interesting to note that there is a correspondence between the string excitation and 
the generalized susceptibilities. For instance the dipole polarizability {dtp) 2 on the world 
line maps to the dilaton excitation of the string. Higher order multipoles map to massive 
string excitations. 



Given the action the two body problem can be solved. We utilize the holomorphic nature 
of the fields to re- write the action as 

5 iBt = ^C B (2H^)(^) 

n 

(18) 

Then using the result (jSJ) we have 

(19) 

The exact result for the partition function follows from treating these operators as mass 
insertions and recognizing that the contractions between two identical world lines lead to 
trivial renormalizations that are pure counter-terms. Thus the only relevant vacuum dia- 
grams are those in which no two mass insertions are adjacent. 

Let us consider the case of two particles. To calculate the free energy we draw all possible 
vacuum diagrams which have no disconnected sub-components. We must always add two 
insertions at a time since for an odd number of insertions there would have to be a self-energy 
graph inserted which we renormalize to zero. Diagrams which are related by a reflection 
symmetry are not distinct. At second order in the interaction the diagram reduces down to 




FIG. 1. All connected vacuum diagrams can be written in term of this topology. The different 
colored dots represent insertions at each of the two particles respectively. 

a simple bubble and the resulting free energy is given by 

oo 2n— 1 - 

~f 3F = Y. II (^^)"-V^y9^J 1 (2i? A ^)7 1 (2E B VO^) 

n=l o=l,3,5.. 
X log(z a - Z a+1 ) \og(z a - Z a+ l) 

(20) 

Each n represents two insertions, one of each type, and z 2n +i = Z\. Upon expanding this 
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expression and differentiating one sets z a = z\ and z a+ i = z 2 . Let us now expand this result 
in R/r to compare with the literature on the lower order contributions. The leading order 
contribution is the dipole-dipole interaction. Expanding out (120]) gives 

- 0F^ d = = R\R 2 B d zl d- zl d Z2 d- Z2 log((*i - z 2 ) log(zi - z 2 )) 



1 



R\Rb (21) 



We also generate non-linear terms such as the case with two dipoles of each type in which 
case i = 2. 

- (3F^ } _ d 2 = R 2 A Rl\ V^A, y/d^I 1 (2R A y/d^)I 1 (2R B y/d^) logfr - z 2 ) \og{ Zl - z 2 ) 
x a/ d Z3 d S:! a/ d Z4 d- ZA h (2R A y/ d Z3 d Z3 ) h (2R B y/ d Z4 d S4 ) log(^ 3 - z A ) \og(z 3 - z A ) 
— f? 4 p 4 

(22) 



All of these perturbative results agree with those in the literature 
read off by going to higher order. 



0,3- 



New results can be 



The case of N holes follows in a similar fashion. If we are interested in the N-body 
potential we have a minimum of N insertions and we may add one vertex insertion at a 
time. 

oo n 

- $F N =J2 E n^ 2 --^ V^^MZRaVdZdZTJ x G(z a - z a+1 , z a - z a+l ) 

n=N V a=l 

(23) 

The second sum is over all possible choice non-equivalent permutation of vertices. Where 
two permutations are in the same equivalence class if they can be rotated into each other. 
Defining G(0) = implies that all choices with two adjacent identical vertices vanish. For 
each choice of vertices, (ii, i 2 , is-..) correspond to the number of vertices of (hole) type 
(1,2,3...) in the particular term in the sum. The expression ( 123]) is easily evaluated to 
arbitrary order in MATHEMATICA, though this particular form is most probably not the 
most efficient since it entails redundancies as compared to the two body results presented 
above. 
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II. CFT SOLUTION FOR TWO HOLES 



The determinant of the Laplacian on Riemann surfaces in many instances can recovered 
by conformally mapping surfaces for which the partition is easily calculated to surfaces 
which present more of a challenge, such as the case of interest in this paper. The techniques 
used here are well known though it seems that the explicit result for the case of interest, 
two spheres in a plane, is absent in the voluminous literature on the subject of functional 
determinants in two dimensions. 



The partition function on the plane with two holes can be calculated by a series of con- 
formal maps starting from the cylinder. One could of course map directly from the cylinder 
to the plane with holes, but it is perhaps simpler to use the annulus as an intermediate step. 
In each step we work in a conformal gauge 



9ab ~ e 2 *g ab , (24) 



in which case the following quantity is an invariant 10|, 



lN = Def N (^) {e ^ IdMKdsy (25) 

Where we have chosen Neumann boundary conditions for the purpose of illustration. L is 
defined as 

L((f>,g)= f dp,K(/>+ [ dsk<t> + \f d(ig a \d a <t>d h <t>), (26) 

J M JdM * J M 

and hatted quantities are defined by the hatted metric, K is the Gaussian curvature on M, 
k is the geodesic curvature on the boundary, and A is the area of M. The prime above the 
determinant denotes the fact that the zero mode is excluded. 



Beginning with the partition function on a cylinder with length / and unit radius. A 
rudimentary calculation yields the result 

1 00 1 

Z = -leT TT L_. (27) 

v 71 = 1 



s 



We then map this result to the annulus using the conformal transformation 



z = r 2 e 



(28) 



with 

2tt/ = ln(r 2 /ri). (29) 
The outer (inner) radii of the annulus are ^(ri). The conformal factor is given by 



-\og{2n zz). 



Given that both spaces are flat we find 

Lann = 7T \og(r 2 /r 1 ) . 



(30) 



(31) 



The annulus partition function is then given by 11] 



z ann = f MM=_ (n/r,)- 1 / 6 



oo 1 
n=l 



(32) 



(33) 



Then we can arrive at case of interest (z plane with two holes) via a conformal map from 
the annulus (w plane). Tha mapping shown in figure ([2]) is given by 
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FIG. 2. Mapping between annulus and twice punctured plane 



w 



z — a 
az — 1 



(34) 
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where 



1 + xix 2 + \/ (x\ - l){xl 
Xi + x 2 



(35) 



R = XiX 2 - 1 - y/(xf - 1)(X 2 2 - 1) ^ 
Xi - X 2 



under the conditions x 2 < a < x%. 

The metric on the annulus is then mapped to the flat z plane 



(a 2 - l) 2 

dzdz = dwdw- r^-. — = -rr (37) 

(aw — iy(aw — \y 



so that 

rh = 

2 \{aw- l) 2 (aw- l) 2 , 

and g = 5 a b is the metric on the flat plane so that K = 0. The annulus is flat (R ~ V 2 </>) so 
IdM K = since its Euler characteristic x — 0. However, the integral 

/» /•27T 

/ dsk(f)= j d6((j>(l,6)-(j>(R,6)) = 2irLog[R) (39) 

does not vanish since is not constant along the boundary. 

The bulk contribution to the anomaly must vanish since the conformal map is an element 
of the SL(2, C) subgroup with vanishing Schwarztian. Hence we are left the the following 
partition function for the plane with two holes. 

- 0F = Log[f[R 2 }} (40) 

where 

oo .. 
n=l v J 

and 

R= l/2(r 2 -r v / r 2 -4-2). (42) 
Expanding in the ratio of the size of the holes relative to the distance between them gives 
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Which agrees with the expansion of (I2U|) . The expression (HU1) certainly has many advantages 
over the result fl20|) . In particular, though the expansion in (1201) has a finite radius of con- 
vergence we can not expand it around to the osculating limit. Whereas the Dedekind func- 
tion which results from (|40[ 41) has multiple representations which can be easily expanded 
around this limit. The result ( 1401) may be used be compared to the so-called "proximity force 



approximation" 
13. 



121 ] which has recently been systematized in terms of a derivative expansion 



III. NON-CONFORMAL THEORIES 

It is interesting to see that the effective theory methodology can still yield exact results 
even for non-conformal kernels. A non-conformally invariant kernel of particular interest is 
the bi-Laplacian which corresponds to a model of a membrane where bending dominates 
tension. In this case there are more operators at each order in R since the equations of 
motion can no longer be used to eliminate operators involving a single Laplacian. In the 
case of the bi-Laplacian, beyond second order in derivatives we have two operators at each 
order in the derivative expansion. At the quadrapole level we have two possible operators. 
A convenient basis to choose is 

2 a = d&WM, 2 b = d z d- z <pd- z d z <p. (44) 

At higher orders we add one power of d Z)Z to each side. So that at n'th order in derivatives 
we have 

Q n+1 = gn+l^gn+1^^ qu+1 = gn g^gn g ^ (45) 

In the case of the bi-Laplacian one must impose two boundary conditions to find unique 
solutions. Which operator has the first non-vanishing susceptibility will depend upon the 
particular choice. The crucial point is that the number of operators at each level in deriva- 



tives does not grow which al 



in the derivative expansion 19]. 



ows again for the matching procedure to carried out to all order 



This fact may at first seem surprising given that the theory is not conformal. However, 
the theory is invariant under the SL(2, C) sub-group of the Virasoro group which is revealed 
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once one makes a field redefinition 15|, |l6] 



(46) 

In addition, higher order kernels can also be solved by the same method. However, as one 
goes to higher orders more (but still finite) operators appear at each order. 

IV. CONCLUSIONS 

Partition functions for on higher genus surfaces can have a rich structure. The boundaries 
introduce new sets of operators which make even free theories highly non-trivial. Models 
which are integrable on simple manifolds, such as cylinders, can be solved on more compli- 
cated surfaces through conformal mappings and the associated anomalous transformations 
of the stress energy in the bulk and boundary q However, mappings to the Euclidean plane 
with multiple boundaries of various shapes can become analytically intractable. 

In this letter we have shown that two dimensional partition functions can also be ex- 
pressed in terms of an infinite set of Wilson coefficients in a world line EFT. For Gaussian 
theories, such as the free bosonic one discussed here, these coefficients are exactly calculable. 
These coefficients which physically correspond to generalized susceptibilities, are equal to 
the weights of string excitations in the coherent state boundary operator in the closed string 
channel. Given these coefficients, the partition function follows from the vacuum energy 
graphs. The result is written as an infinite sum which is clearly inferior to the infinite prod- 
uct form which arise in the context of CFT methodology. For the case of two spherical holes 
we have shown the equality of the CFT and EFT results expanding to arbitrary order in 
the separation between the holes. 

However, as the number of defects increases the only impediment to solving the theory 
using the world-line methodology is that the series involved grows more cumbersome. Here 
we only presented the sums in their most raw form, and it is surely possible to simplify 



them further 



20| | . Furthermore, it is also possible to determine the Wilson coefficients for 



deformed holes using the fact that can take advantage of rudimentary conformal transforma- 



tion techniques to calculate the infinite set of Wilson coefficients for ellipses [171] • We have 



The partition function for sections of spheres have been calculated in this fashion in 
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also shown that it is also possible to solve for the Wilson coefficients for the non-conformal 
bi-Laplacian kernel. 

Finally it is interesting to ask how these EFT techniques can be applied to interacting 
theories. For any conformal theory, the conformal mapping technique is still applicable, and 
here is no reason that one could not find the fluctuation induced force say for an Ising model 



181 ] . However, again if one is interested in more complicated geometries (such as multiple 



holes) the mappings will become analytically intractable. However, given that distortions of 
single holes to other shapes, such as ellipses, is fairly simple, the matching procedure in the 
effective theory technique is still relatively straightforward. Thus, given that for a simple 
geometry one can solve for the Wilson coefficient, at least for an integrable model, one can 
also determine the coefficients for more complicated shapes. The free energy for multiple 
holes can then be calculated via the vacuum energy diagrams as shown above, at least for 
theories which are integrable. 
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